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Abstract 

This  paper  presents  an  economic  model  of  a  PEM  fuel  cell  power  plant  (FCPP).  The  model  includes  the  operational  cost,  thermal  recovery, 
power  trade  with  the  local  grid,  and  hydrogen  production.  The  model  is  used  to  determine  the  optimal  operational  strategy,  which  yields  the 
minimum  operating  cost.  The  optimal  operational  strategy  is  achieved  through  estimation  of  the  following:  hourly  generated  power,  thermal 
power  recovered  from  the  FCPP,  power  trade  with  the  local  grid,  and  hydrogen  production.  An  evolutionary  programming-based  technique 
is  used  to  solve  for  the  optimal  operational  strategy.  The  model  is  tested  using  different  seasonal  load  demands.  The  results  illustrate  the 
impact  of  producing  hydrogen  on  the  operational  strategies  of  the  FCPP  when  subjected  to  seasonal  load  variation.  Results  are  encouraging 
and  indicate  viability  of  the  proposed  model. 
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1.  Introduction 

In  recent  years  distributed  sources  have  received  consider¬ 
able  attention  because  of  the  improvements  in  power  quality, 
reliability,  portability,  and  environmental  emissions.  Among 
the  various  types  of  distributed  generators,  fuel  cell  power 
plants  (FCPPs)  have  been  the  focus  of  interest  since  such 
plants  are  capable  of  producing  electricity,  heat,  and  hydro¬ 
gen.  Due  to  low  working  temperature  (80-100  °C),  fast  start 
up,  extremely  low  emission,  and  very  low  noise  PEM  FCPPs 
are  the  best  candidates  for  residential  and  isolated  load  appli¬ 
cations.  With  cost  effective  operational  strategies,  the  use  of 
FCPPs  is  expected  to  become  widespread  in  the  near  future, 
in  spite  of  their  current  high  capital  cost. 

Fuel  cell  economics  and  economical  aspects  have  been 
presented  in  the  literature  [1-5].  In  [1,2]  an  economic  model 
has  been  introduced  to  estimate  the  optimal  output  power 
from  the  FCPP  while  satisfying  system  operational  con¬ 
straints.  This  simple  model  considers  only  the  possibility  of 
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selling  and  buying  energy  from  the  local  grid,  and  the  utiliza¬ 
tion  of  thermal  power  output  from  the  FCPP. 

In  this  paper  the  model  in  [1,2]  has  been  extended  to 
include  hydrogen  production  of  the  FCPP.  The  hydrogen  pro¬ 
duction  is  considered  at  times  when  the  electrical  load  is 
lower  than  the  maximum  capacity  of  the  FCPP.  Introduction 
of  hydrogen  generation  in  the  model  brings  to  completion  all 
the  operational  aspects  of  the  FCPP  that  have  an  impact  on 
the  reduction  of  the  overall  system  cost  and  increase  in  the 
FCPP  efficiency. 

The  economic  model  is  represented  as  a  cost  optimization 
problem  subject  to  system  and  operational  constraints.  To 
estimate  the  daily  optimal  operational  strategy  for  the  FCPP 
a  hybrid  technique  based  on  evolutionary  programming  (EP) 
and  Hill-Climbing  (HC)  method  [1,6]  is  used.  The  evolution¬ 
ary  programming  is  employed  to  search  for  the  near  optimal 
solution  while  the  HC  method  is  used  to  ensure  feasibility 
during  the  solution  process. 

The  paper  is  organized  as  follows:  Section  2  introduces  an 
economic  model  for  an  FCPP  system.  Section  3  presents  the 
solution  methodology.  Test  results  are  presented  in  Section  4 
and  Section  5  presents  the  conclusions. 
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2.  Fuel  cell  system  economic  model  development 

In  this  model  many  different  strategies  are  developed  to 
handle  excess  electrical  and  thermal  energy  and  hydrogen 
production. 

2.7.  Fuel  cell  economic  model 

In  [1,2],  the  authors  introduced  a  mathematical  model  for 
the  FCPP  operational  cost.  In  this  paper,  the  model  has  been 
extended  to  include  the  economic  aspects  of  hydrogen  gener¬ 
ation.  The  model  considers  electrical  power  generation,  ther¬ 
mal  power  recovery,  and  hydrogen  production.  The  model 
given  below,  represents  the  operational  cost  in  US$day-1, 
and  can  be  summarized  as  follows: 


MUT:  minimum  up-time  (number  of  intervals); 

MDT:  minimum  down-time  (number  of  intervals); 

U :  FCPP  on-off  status,  U  =  1  for  running,  U  =  0  for  stop¬ 

ping; 

A^imx;  maximum  number  of  start-stop  events; 

^start-stop •  number  of  start- stop  events; 

OM:  operation  and  maintenance  cost; 

CHs:  hydrogen  selling  price  (US$  kg-1); 

Puf  the  equivalent  electric  power  for  hydrogen  produc¬ 
tion  (kW); 

F:  a  conversion  factor  (kg  of  hydrogen/kW  of  electric 

power),  where  F=  1.05  x  10-8/nceii  and  i>cen  is  cell 
operating  voltage,  uceii  =  0.6  volt; 

Pjj:  total  power  produced  at  interval  j,  where  Pjj  = 

Pj  +  Pa  +  PUj 


min 


Pj  +  +  A  I;  \ 

1]  ) 


+  CeijPT  y^  maxdei.  /  -  Pj,  0) 

j 


Ca,sT  ^  ma x(Pj  -  Le \j,  0) 
j 


+  CnlT  yy  max(Lth,7-  -  Ah,;,  0)  +  a  +  p(l  -  e  r°ff/r)  +  OM 
j 


Ch.tY^PbjF 

j 


Subject  to: 

pmin  <  pT  <  pmnx 

(2) 

Pj  -  Pj- 1  <  aa 

(3) 

pj- 1  “  pi  A  A  Ad 

(4) 

(7™,  -  MUTXt/,-1  -  Uj)  >  0.0 

(5) 

(TfJ  -  MDT )(Uj  -  Uj- 1)  >  0.0 

(6) 

„  ^  \  j  max 

n  start- stop  _  l  s 

(7) 

where 

Cni :  price  of  natural  gas  for  FCPP  (US$  kW  h-1); 

T:  length  of  time  interval  (h); 

Pj  :  electrical  power  produced  at  interval  j  (kW)  less  the 

power  for  auxiliary  devices; 

Pa:  power  for  auxiliary  devices  (kW); 

iy:  fuel  cell  efficiency  at  interval  j ; 

Cei?p :  tariff  for  purchasing  electricity  (US$  kW  h- 1 ); 

Cels :  tariff  for  selling  electricity  (US $  kW  h- 1 ) ; 

Le\j :  electrical  load  demand  at  interval  j  (kW); 

Cn 2 :  fuel  price  for  residential  loads  (US$  kW  h- 1 ); 

Lthj:  thermal  load  demand  at  interval  j  (kW); 

Pthj  •  thermal  load  produced  at  interval  j  (kW); 
a,  f3 :  hot  and  cold  start  up  cost,  respectively; 

toft"  time  the  FCPP  has  been  off  (h); 
r:  fuel  cell  cooling  time  constant  (h); 

Pmin:  minimum  limit  of  generating  power  (kW); 

pmax.  maximum  limit  of  generating  power  (kW); 

A Pu:  upper  limit  of  the  ramp  rate; 

A Pj)  \  lower  limit  of  the  ramp  rate; 

7™:  FCPP  on-time  (number  of  intervals); 

T°ff:  FCPP  off-time  (number  of  intervals); 


First  term  of  the  objective  function  is  the  daily  overall  fuel 
cost  for  the  FCPP  (US$).  Second  term  is  the  daily  cost  of  elec¬ 
trical  energy  purchased  if  the  demand  exceeds  the  electrical 
energy  produced  (US$).  Third  term  is  the  daily  income  from 
the  electrical  energy  sold  if  the  electrical  energy  produced 
exceeds  the  demand  (US$).  The  fourth  term  is  the  daily  cost 
of  purchased  gas  for  residential  thermal  loads  if  the  thermal 
energy  produced  is  not  enough  to  meet  the  thermal  energy 
demand  (US$).  The  fifth  term  is  the  start  up  cost  (US$).  The 
sixth  term  is  the  operation  and  maintenance  cost  of  the  FCPP 
(US$).  The  last  term  is  the  daily  income  from  the  sale  of 
hydrogen  generated  by  the  FCPP  (US$). 

2.2.  Excess  electrical  energy  strategy 

In  the  model,  excess  electrical  energy  can  be  sold  to  the 
local  grid  according  to  system  economics.  Under  certain  con¬ 
ditions,  it  might  be  beneficial  to  the  system,  in  terms  of  the 
overall  cost,  to  buy  energy  from  the  local  grid  to  satisfy  the 
load  requirements  while  using  the  unused  capacity  to  produce 
hydrogen  and  thermal  energy. 

2.3.  Hydrogen  production  strategy 

The  hydrogen  production  strategy  is  based  on  the  differ¬ 
ence  between  the  maximum  capacity  of  the  FCPP  and  the 
generated  electric  power  at  each  interval. 

To  include  the  hydrogen  in  the  FCPP  model,  an  equiva¬ 
lent  electric  power  for  the  generated  hydrogen  at  each  interval 
is  considered  Pnj.  The  equivalent  electric  power  is  consid¬ 
ered  at  the  fuel  cell  stack  output  as  shown  in  Fig.  1.  Fig.  1 
(upper  figure)  reflects  the  electric  power  output  and  hydro¬ 
gen  output  locations  in  the  FCPP  stages.  Fig.  1  (lower  figure) 
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Fig.  1.  Hydrogen  insertion  in  the  FCPP  model. 

shows  the  location  of  Puj  in  the  FCPP.  Considering  Puj  at 
the  stack  terminals  makes  it  possible  to  quantify  the  produc¬ 
tion  of  hydrogen  (kgs-1)  in  terms  of  electrical  energy.  The 
hydrogen  production  in  kgs-1  can  be  calculated  using  Piy 
as  follows  [7] : 

(H2) amount  =  1.05  x  10-8— ^  (8) 

Fcell 

In  this  paper  two  strategies  for  generating  hydrogen  are 
adopted: 

2.3.1.  Strategy  1 

The  amount  of  the  generated  hydrogen  is  equal  to  the  dif¬ 
ference  between  the  maximum  capacity  and  the  generated 
power  level  of  the  FCPP  as  given  in  Eq.  (9). 

Ph  j  —  Pmax  Pj  P a  (9) 

This  strategy  ensures  rated  capacity  operation  of  the 
reformer  of  the  FCPP  all  the  time.  Rated  capacity  output 
allows  the  reformer  to  be  more  stable  and  work  at  high  effi¬ 
ciency. 


Fig.  2.  Performance  curves  of  the  FCPP. 


energy  approximately  equal  to  the  electrical  energy  [8].  In 
[8],  efficiency  and  thermal  energy  to  electrical  energy  ratio 
curves  have  been  developed  (Fig.  2).  These  curves  approx¬ 
imate  the  efficiency  and  the  thermal  output  of  the  FCPP. 
The  efficiency  and  the  thermal  energy  to  electrical  energy 
ratio  are  functions  of  the  part  load  ratio  (equal  to  electrical 
generated  power/maximum  power).  Mathematical  expres¬ 
sions  to  approximate  the  curves  have  been  developed  in  [8] 
as  follows: 

For  PLRy  <  0.05 

tj  j  =  0.2716,  rTEj  =  0.6801  (10) 

For  PLRy  >0.05 

y]j  =  0.9033PLR]  -  2.9996PLR]  +  3.6503PLR? 

-  2.0704PLRy  +  0.4623PLR,  +  0.3747  (11) 

rTEJ  =  1.0785PLR]  -  1 .9739PLR;]  +  1.5005PLRj 

— 0.2817PLR,  +  0.6838  (12) 

where  ry  is  the  FCPP  efficiency,  PLR  the  part  load  ratio,  rEE 
the  thermal  energy  to  electrical  energy  ratio. 

The  thermal  power  recovered  from  the  fuel  cell  according 
to  electric  and  hydrogen  power  outputs  can  be  calculated  as 
follows: 


2.3.2.  Strategy  2 

In  this  strategy  hydrogen  production  can  vary  between 
zero  and  the  difference  between  the  maximum  capacity  and 
the  generated  electric  power.  This  strategy  allows  the  model 
to  produce  hydrogen  when  it  is  more  profitable  to  the  overall 
cost  of  FCPP  operation. 

2.4.  Recovered  thermal  energy  strategy 

The  FCPP  operates  with  approximately  36%  efficiency. 
The  efficiency  is  slightly  higher  at  low  load  compared  to 
full  load  operation.  At  full  load,  the  FCPP  produces  thermal 


Ah,  j  =  7Te(Phj  +  Pj  +  Pa) 


(13) 


The  efficiency  of  the  FCPP  given  in  Fig.  2  is  based  on  the 
electrical  output  power  versus  the  input  gas  power.  Neglect¬ 
ing  the  thermal  power  in  efficiency  calculation  results  in 
efficiency  range  of  30-40%  approximately.  On  the  other 
hand,  including  the  utilized  thermal  power  recovered  from 
the  FCPP  and  hydrogen  production  enhances  the  overall  effi¬ 
ciency  considerably.  The  overall  efficiency  can  be  calculated 
as  follows: 


7?  overall  j  — 


Pj  +  Pa  +  Puj  +  min(Lthj,  Pthj) 
t(P j  +  Pa  +  Ph  j)  /  0  j] 


(14) 
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3.  Evolutionary  programming  (EP)-based  solution 
methodology 

Evolutionary  programming  can  be  traced  back  to  the 
early  1950s  when  Turing  discovered  a  relationship  between 
machine  learning  and  evolution  [9-11].  Later,  Bremermann, 
Box,  Friedberg,  and  others  developed  evolutionary  compu¬ 
tation  as  a  tool  for  machine  learning  and  optimization.  Great 
attention  was  given  to  EP  as  a  powerful  tool  when  Fogal,  Bur- 
gin,  Atmar,  and  others  used  it  to  predict  the  events  of  finite 
state  machines  on  the  basis  of  old  observations.  During  the 
1980s  evolutionary  programming,  with  advances  in  computer 
technology,  was  used  to  solve  difficult  real-world  optimiza¬ 
tion  problems.  In  the  power  systems  area,  EP  has  been  used 
to  solve  a  number  of  power  systems  problems  [11]. 

Evolutionary  programming  is  a  search  optimization 
method.  It  moves  from  one  solution  to  another  using  a  prob¬ 
abilistic  search  technique.  Evolutionary  programming  starts 
with  random  individuals.  Each  individual  represents  a  com¬ 
plete  solution  for  the  problem  under  study.  The  individuals 
are  moved  from  one  generation  (or  iteration)  to  the  other  after 
passing  through  two  main  steps,  mutation  and  competition. 
During  a  mutation  step  a  new  individual  is  produced  when  a 
Gaussian  random  variable  with  uniform  probability  is  added 
to  the  current  individual.  The  competition  step  is  a  probabilis¬ 
tic  selection  scheme  used  to  assign  a  weight  to  each  individual 
according  to  a  comparison  between  current  individual  and  a 
randomly  chosen  one.  It  may  happen  that  the  new  solution 
is  infeasible.  Therefore,  using  EP  alone  may  require  a  long 
time  to  reach  the  optimal  solution  or  it  may  get  trapped  in 
a  local  optimum.  This  limitation  was  overcome  by  the  use 
of  the  HC  technique  [12]  to  move  new  infeasible  solutions 
into  the  feasible  region.  The  following  algorithm  details  the 
proposed  approach  to  solve  the  problem: 

1 .  Generate  initial  random  solutions  for  the  output  power 
from  the  FCPP  at  each  interval. 

Si  =  {v},  i  =  1, . . . ,  m  (15) 

where  v  is  a  set  of  output  power  from  the  FCPP  at  each 
interval;  m  the  number  of  individual  in  the  current  gen¬ 
eration. 

The  random  solution  is  expected  to  satisfy  the  system 
constraints. 

2.  For  each  individual  in  the  current  generation,  calculate 
the  objective  function  value  using  (1). 

3.  Mutate  each  individual  and  assign  it  to  Si+m  according 
to  (16). 

—  Si  +  N(0,Piv(Si)  +  Zi)  (16) 

where  Si  is  zth  individual;  k  the  number  of  generating 
units  to  be  maintained  in  the  current  individual;  A(/x,  a1 2 3) 
the  Gaussian  random  variable  with  mean  p  and  variance 
cr2;  Pi  a  constant  to  scale  n(S));  Zj  an  offset  to  guarantee 
a  minimum  amount  of  variance. 


4.  Check  the  feasibility  of  each  new  individual  against  the 
constraints.  If  there  is  no  violation  go  to  step  5.  Otherwise 
go  to  step  6. 

5.  Calculate  the  objective  function  value  for  the  feasible 
solution  using  (1)  and  go  to  step  7. 

6.  Use  the  Hill-Climbing  algorithm  to  drive  the  infeasible 
individuals  into  feasibility.  If  no  feasible  solution  can  be 
found  go  to  step  3. 

7.  Assign  a  fitness  score  v(Si)  to  each  individual  Si+m  (i=  1 , 
. . .,  2m).  The  score  is  assigned  equal  to  the  cost  function. 

8.  Using  Eq.  (17),  calculate  a  weight  Wj  for  each  individ¬ 
ual  Si,  i=  1,  . . .,  2m.  These  weights  are  to  be  calculated 
during  a  random  competition  between  individuals  based 
on  the  objective  function  value. 

N 

Wi  =  ]T  Wij  (17) 

7=1 

where  A  is  a  randomly  generated  competition  number; 
Wq  either  0  or  1  depending  on  the  competition  of  the 
individual  with  another  individual  selected  randomly 
from  the  population.  The  value  of  Wq  can  be  calculated 
as  follows: 

1  if  v(Si)<v(Sp) 

(to) 

0  otherwise 

where 

p  =  [2mu\  +  1],  p  ^  i  and  u\  ~  U( 0,  1) 

9.  Rank  the  solution  S((i  =  1,  . . .,  2m)  in  descending  order 
according  to  their  values  of  Wi  (if  more  than  one  solution 
has  the  same  W,  use  the  actual  score  of  v(Si)  to  rank 
them).  Use  the  first  m  solutions  along  with  their  score 
values  v(Si)  as  a  new  generation  for  the  potential  optimal 
solution. 

10.  Check  for  convergence.  Criteria  used  for  convergence 
include  the  maximum  generation  number  and  the  aver¬ 
age/maximum  fitness  ratio  being  less  than  a  predeter¬ 
mined  small  value.  If  convergence  is  achieved,  stop; 
otherwise  go  to  step  3. 

4.  Tests  and  results 

The  proposed  model  has  been  applied  to  a  250  kW  grid- 
parallel  FCPP  that  supplies  a  residential  neighborhood.  The 
IEEE-RTS  load  profile  with  a  peak  of  250  kW  [13]  is  used  to 
simulate  the  hourly  electrical  load  profile  of  the  system.  In 
this  test  system,  the  weekly,  daily  and  hourly  peak  load  values 
are  given  in  percent  of  annual,  weekly  and  daily  peak  loads, 
respectively.  The  winter  hot  water  usage  and  space  heating 
load  for  Atlanta,  Georgia  [8]  is  considered  to  represent  the 
thermal  load  profile.  Due  to  the  lack  of  thermal  load  infor¬ 
mation  for  the  summer  and  spring/fall,  thermal  load  data  are 
estimated  from  the  available  winter  data.  The  thermal  load 
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Table  1 

FCPP  and  evolutionary  program  parameters 

Maximum  limit  of  generating  power,  Pmax  (kW)  250 

Minimum  limit  of  generating  power,  Pmm  (kW)  0.0 

Length  of  time  interval,  T  (h)  0.25 

Upper  limit  of  the  ramp  rate,  A Pu  (kW  s-1)  20 

Lower  limit  of  the  ramp  rate,  APd  (kW  s-1)  25 

Price  of  natural  gas  for  FCPP,  Cn\  (US$  kW  h-  1 )  0.04 

Tariff  for  purchasing  electricity,  Ce i,p  (US$  kW  h- 1 )  0.13 

Tariff  for  selling  electricity,  Cei,s  (US$  kWh-1)  0.08 

Fuel  price  for  residential  loads,  Cn 2  (US$  kW  h-1 )  0.06 

Hydrogen  selling  price,  Chs  (US$kg-1)  1.80 

Hot  start  up  cost,  a  (US$)  0.05 

Cold  start  up  cost,  f3  (US$)  0.15 

The  fuel  cell  cooling  time  constant,  z  (h)  0.75 

Minimum  up-time,  MUT  (number  of  intervals)  2 

Minimum  down-time,  MDT  (number  of  intervals)  2 

Maximum  number  of  start-stop  time,  /Vmax  5 

Maximum  number  of  evolutionary  generation  20000 

Number  of  individuals  150 


is  used  along  with  the  electrical  load  profile  to  simulate  total 
hourly  operation  of  the  FCPP.  In  the  following  case  stud¬ 
ies,  the  optimum  operational  cost  is  evaluated  and  compared 
with  a  Base  Case.  The  gas  prices,  hydrogen  selling  price,  and 
FCPP/EP  parameters  for  all  test  cases  are  given  in  Table  1. 

4.1.  Base  Case 

In  this  case  the  above  thermal  and  electrical  loads  are 
used  to  estimate  the  optimum  operational  strategy  for  FCPP 
operation  without  hydrogen  production.  The  obtained  cost 
components  for  different  seasons  are  given  in  Table  2. 

4.2.  Case  1 

In  Case  1,  Strategy  1  is  tested  with  different  seasonal  ther¬ 
mal  and  electrical  loads.  The  results  for  the  cost  components 
for  different  seasons  are  given  in  Table  3.  It  is  evident  from 
comparison  of  Tables  2  and  3  that  hydrogen  production  saves 
US$  62.28,  US$  27.42  and  US$  15.15  daily  in  the  winter, 
summer  and  spring/fall  cases,  respectively.  The  total  sav¬ 
ing  for  the  winter,  summer  and  spring/fall  seasons  are  US  $ 
5698.62,  US$  2508.93  and  US$  2772.45,  respectively.  Ther¬ 
mal  load  impacts  hydrogen  production  significantly  as  shown 
in  the  spring/fall  case.  Fig.  3  shows  electrical  power  genera¬ 
tion  and  load  for  the  spring/fall  case.  In  this  case,  there  was  no 
electrical  power  trade  with  the  local  network.  Fig.  4  shows 

Table  2 


Cost  component  for  Base  Case 


Cost  components  (US$) 

Winter 

Summer 

Spring/fall 

Daily  fuel  cost 

613.91 

602.90 

539.72 

Daily  cost  of  purchased  electricity 

0.00 

0.00 

2.31 

Daily  profit  from  electricity  sold 

76.60 

115.82 

99.48 

Daily  cost  of  residential  natural  gas 

149.79 

29.16 

45.55 

Operation  and  maintenance  cost 

26.37 

26.04 

23.79 

Total  cost 

713.47 

542.27 

511.88 

Table  3 


Cost  component  for  Strategy  1 


Cost  components  (US$) 

Winter 

Summer 

Spring/fall 

Daily  fuel  cost 

748.60 

748.60 

748.60 

Daily  cost  of  purchased  electricity 

0.00 

2.38 

0.00 

Daily  profit  from  electricity  sold 

0.00 

0.00 

0.00 

Daily  cost  of  residential  natural  gas 

71.87 

1.24 

10.39 

Operation  and  maintenance  cost 

21.58 

18.71 

17.66 

Daily  hydrogen  selling  profit 

190.86 

256.07 

279.92 

Total  cost 

651.19 

514.85 

496.73 

hydrogen  production,  and  thermal  power  load  and  genera¬ 
tion  for  the  spring/fall  case.  Examining  the  thermal  power 
load  in  Fig.  4b  shows  that  the  system  experiences  a  low  ther¬ 
mal  load  between  8:15  a.m.  and  6:00  p.m.  During  this  period 
the  FCPP  is  seen  to  produce  more  thermal  power  than  the 
thermal  load  demand.  This  is  because  the  FCPP  is  forced  to 
produce  hydrogen  equal  to  the  difference  between  the  rated 
capacity  and  the  generated  electric  power.  Producing  ther¬ 
mal  power  more  than  the  requirements  of  the  thermal  load 
makes  the  system  lose  money.  Comparing  the  total  cost  for 
spring/fall  from  Tables  2  and  3  shows  that  the  system  only 
saves  US$  15.15  daily  when  producing  hydrogen  during  low 
thermal  demand  periods. 

4.3.  Case  2 

In  this  case  hydrogen  production  is  based  on  the  sys¬ 
tem  economics  as  in  Strategy  2.  Hydrogen  production  varies 
between  zero  and  the  difference  between  the  rated  capacity 
of  the  FCPP  and  the  generated  electric  power. 

Using  the  thermal  and  electrical  load  profiles  for  differ¬ 
ent  seasons,  the  model  gives  the  optimum  cost  as  shown  in 
Table  4.  Figs.  5-10  show  the  following:  electrical  power  trade 
with  the  grid;  electrical  power  load  and  generation;  hydrogen 
production;  and  thermal  load  and  generation  for  different  sea¬ 
sons.  In  the  winter  case,  there  was  no  electrical  power  trade 
with  the  local  network.  Comparing  different  cost  components 
in  Tables  3  and  4  shows  that  Strategies  1  and  2  yield  the  same 
daily  cost  for  the  winter  season  because  of  the  high  ther¬ 
mal  demand.  For  summer  and  spring/fall  cases  Strategy  2 
gives  lower  operating  cost.  In  the  summer  season,  Strategy 
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Fig.  3.  Spring/fall  electrical  load  and  power  generation. 
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Fig.  4.  (a)  Spring/fall  hydrogen  production;  (b)  spring/fall  thermal  load  and  generation. 
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Table  4 


Cost  component  for  Strategy  2 


Cost  components  (US$) 

Winter 

Summer 

Spring/fall 

Daily  fuel  cost 

748.60 

652.08 

596.07 

Daily  cost  of  purchased  electricity 

0.00 

32.99 

40.92 

Daily  profit  from  electricity  sold 

0.00 

12.32 

4.08 

Daily  cost  of  residential  natural  gas 

71.87 

1.33 

10.46 

Operation  and  maintenance  cost 

21.58 

18.30 

16.34 

Daily  hydrogen  selling  profit 

190.86 

206.79 

202.27 

Total  cost 

651.19 

485.59 

457.44 

Fig.  5.  Winter  electrical  load  and  power  generation. 


CM  CM 


(a) 


Time  (hh:mm) 


(b) 


Time  (hh:mm) 


Fig.  6.  (a)  Winter  hydrogen  production;  (b)  winter  thermal  load  and  generation. 
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Fig.  7.  (a)  Summer  electrical  power  trade  with  the  grid;  (b)  summer  electrical  load  and  power  generation. 
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Fig.  8.  (a)  Summer  hydrogen  production;  (b)  summer  thermal  load  and  generation. 
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Fig.  9.  (a)  Spring/fall  electrical  power  trade  with  the  grid;  (b)  spring/fall  electrical  load  and  power  generation. 


2  saves  US$  29.26  daily  compared  to  Strategy  1,  and  US$ 
56.68  compared  to  the  Base  Case  (no  hydrogen  production); 
which  totals  to  US$  2677.29  and  US$  5186.22  for  the  sum¬ 
mer  season.  Lower  cost  in  Strategy  2  is  because  of  the  fact 
that  Strategy  2  produces  enough  hydrogen  that  is  profitable  to 
the  system  while  satisfying  the  thermal  power  requirements. 
This  can  be  seen  from  Fig.  8b  where  the  thermal  load  and 
generation  are  almost  identical  except  for  the  peak  values. 

In  spring/fall  season,  Strategy  2  saves  US$  39.29  daily 
compared  to  Strategy  1  and  US$  54.44  daily  compared  to  the 
Base  Case,  which  sum  to  US$  7190.07  and  US$  9962.52  per 
spring/fall  season,  respectively.  In  spring/fall  season,  Strat¬ 


(a)  Time  (hh:mm) 


egy  2  buys  power  from  the  grid  during  the  low  thermal  power 
demand  period  (8:15  a.m.-6:00  p.m.).  In  this  period  the  FCPP 
produces  low  electrical  power  and  hydrogen  so  as  to  produce 
enough  thermal  power  to  satisfy  the  thermal  load  require¬ 
ments.  This  can  be  seen  in  Figs.  9a, b,  and  10a,b. 

In  summary,  the  daily  cost  for  different  seasons  and  strate¬ 
gies  are  shown  in  Table  5.  Comparing  costs  for  Strategies  1 
and  2  reveals  that  Strategy  2  gives  a  lower  cost  over  the  entire 
year.  Using  Eq.  (14),  the  overall  efficiency  for  the  Base  Case, 
and  Strategies  1  and  2  for  the  spring/fall  season  are  given 
in  Fig.  11.  Also  included  in  the  figure  is  the  basic  efficiency 
curve  (without  considering  hydrogen  and  the  recovered  ther- 
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Fig.  10.  (a)  Spring/fall  hydrogen  production;  (b)  spring/fall  thermal  load  and  generation. 
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Table  5 


Cost  summary  (US$) 


Base  Case 

Strategy  1 

Strategy  2 

Cost 

Saving 

Cost 

Saving 

Cost 

Saving 

Winter 

713.47 

0.0 

651.19 

62.28 

651.19 

62.28 

Summer 

542.27 

0.0 

514.85 

27.42 

485.59 

56.68 

Spring/fall 

511.88 

0.0 

496.73 

15.15 

457.44 

54.44 

Time  (hh:mm) 


Fig.  11.  FCPP  efficiency  plots  for  spring/fall  season. 

mal  energy)  for  the  same  season.  It  is  clear  that  considering 
the  hydrogen  and  the  recovered  thermal  energy  enhances 
the  overall  efficiency  as  shown  in  Fig.  11.  Comparison  of 
the  different  efficiency  values  shows  that  Strategy  1  overall 
efficiency  is  lower  than  the  basic  efficiency  value  when  the 
system  experiences  low  thermal  power  demand.  Strategy  1 
efficiency  is  lower  than  that  of  the  Base  Case  and  Strategy 
2  efficiency  at  all  load  conditions.  Strategy  2  gives  lower 
efficiency  than  the  Base  Case  at  high  thermal  load  demand 
intervals.  At  low  thermal  demand  periods  Strategy  2  effi¬ 
ciency  is  higher  than  the  Base  Case.  Although  the  efficiency 
of  Strategy  2  is  lower  than  that  of  Base  Case,  Strategy  2  yields 
a  lower  operating  cost  as  described  previously. 

5.  Conclusions 

In  this  paper,  the  impact  of  hydrogen  production  on  the 
optimal  cost  of  operation  of  a  PEM  FCPP  operating  in  a 
grid-parallel  mode  is  presented.  The  economic  model  of  the 
operational  cost  of  the  FCPP  has  been  developed  which 
includes  power  trade  with  the  local  grid,  thermal  recovery 
and  hydrogen  production.  The  paper  offers  practical  con¬ 
cepts  concerning  operational  cost  modeling  of  the  FCPP.  The 
model  incorporates  two  different  strategies  for  hydrogen  pro¬ 
duction.  The  two  strategies  were  evaluated  using  IEEE  test 
system  load  profiles  for  different  seasons.  From  the  results, 
it  can  be  concluded  that  while  Strategy  2  offers  lower  overall 
efficiency  of  the  FCPP  than  the  Base  Case,  Strategy  2  yields 
lower  operating  cost. 


The  figures  presented  in  this  paper  are  based  on  generic 
load  profiles.  Therefore,  region-specific  load  profiles  would 
yield  results  that  necessarily  differ  from  those  presented  in 
this  paper. 

Acknowledgment 

This  research  was  supported  by  a  grant  from  the  Depart¬ 
ment  of  Energy  (DE-FG02-02ER63376). 


References 

[1]  M.Y.  El-Sharkh,  A.  Rahman,  M.S.  Alam,  Evolutionary 
programming-based  methodology  for  economical  output  power 
from  PEM  fuel  cell  for  micro-grid  application,  J.  Power  Sources 
139  (1-2)  (2005)  165-169. 

[2]  S.  Ahmed,  M.  Azmy,  I.  Erlich,  Online  optimal  management  of  PEM 
fuel  cells  using  neural  networks,  IEEE  Trans.  Power  Delivery  20  (2) 
(2005)  1051-1058. 

[3]  F.  Barbir,  T.  Gomez,  Efficiency  and  economics  of  proton  exchange 
membrane  (PEM)  fuel  cell,  Int.  J.  Hydrogen  Energy  21  (10)  (1996) 
891-901. 

[4]  Y.  Zoka,  H.  Sasaki,  J.  Kubokawa,  R.  Yokoyama,  H.  Tanaka,  An 
optimal  deployment  of  fuel  cells  in  distribution  systems  by  using 
genetic  algorithms,  IEEE  Int.  Conf.  Evolut.  Comput.  1  (1996),  479, 
484. 

[5]  G.  Erdmann,  Future  economics  of  the  fuel  cell  housing  market,  Int. 
J.  Hydrogen  Energy  28  (2003)  685-694. 

[6]  M.Y.  El-Sharkh,  A.A.  El-Keib,  Maintenance  scheduling  of  power 
system  generation  and  transmission  using  fuzzy  evolutionary  pro¬ 
gramming,  IEEE  Trans.  Power  Syst.  18  (2)  (2003). 

[7]  J.E.  Larmine,  A.  Dicks,  Fuel  Cell  Systems  Explained,  second  ed., 
John  Wiley  and  Sons,  2000. 

[8]  M.B.  Gunes,  Investigation  of  a  fuel  cell  based  total  energy  system 
for  residential  applications,  Master  of  Science  thesis,  Department 
of  Mechanical  Engineering,  Virginia  Polytechnic  Institute  and  State 
University,  2001. 

[9]  D.B.  Fogel,  Evolutionary  Computation  toward  a  New  Philosophy  of 
Machine  Intelligence,  second  ed.,  Wiley-IEEE  Press,  1999. 

[10]  T.  Back,  U.  Hammel,  H.P.  Schwefel,  Evolutionary  computation:  com¬ 
ments  on  the  history  and  current  state,  IEEE  Trans.  Evolut.  Comput. 
1  (1)  (1997)  3-17. 

[11]  V.  Miranda,  D.  Srinivasan,  L.M.  Proenca,  Evolutionary  computation 
in  power  system,  Electrical  Power  Energy  Syst.  20  (2)  (1998)  89-98. 

[12]  PH.  Winston,  Artificial  Intelligence,  third  ed.,  Addison- Wesley  Pub¬ 
lishing  Company,  1993. 

[13]  Reliability  Test  System  Task  Force  of  the  Application  of  Probability 
Methods  Subcommittee,  IEEE  Reliability  Test  System,  IEEE  Trans. 
Power  Apparatus  Syst.  PAS-98  (6)  (1979)  2047-2054. 

M.Y.  El-Sharkh  received  his  BSc  and  MSc  from  Ain  Shams  University, 
Cairo,  Egypt  in  1987,  and  1990,  respectively,  all  in  electrical  engineer¬ 
ing.  He  obtained  the  PhD  degree  from  the  University  of  Alabama  at 
Tuscaloosa  in  2002.  His  research  interests  include  using  the  new  opti¬ 
mization  techniques  in  solving  the  power  system  operation  problems, 
power  system  stability  and  alternate  sources  of  energy. 

M.  Tanrioven  received  his  BSc,  MSc  and  PhD  from  Yildiz  Technical 
University,  Istanbul,  Turkey  in  1993,  1996,  and  2000,  respectively,  all  in 
electrical  engineering.  His  research  interests  include  power  system,  power 
quality  and  harmonics,  fuzzy  systems,  alternate  sources  of  energy.  He  is 
currently  engaged  in  post-doctoral  research  at  the  University  of  South 
Alabama,  Mobile,  USA. 


144 


M.Y.  El-Sharkh  et  al.  /  Journal  of  Power  Sources  153  (2006)  136-144 


A.  Rahman  received  the  PhD  degree  in  electrical  engineering  from  the 
University  of  Kentucky  in  1974.  He  is  currently  an  associate  professor  of 
electrical  and  computer  engineering  at  the  University  of  South  Alabama. 
His  research  interests  are  in  the  areas  of  electrical  machines,  power  elec¬ 
tronics  and  power  systems.  He  is  a  life  senior  member  of  IEEE  and 
member  of  ASEE. 


M.  Alam  is  a  professor  of  electrical  and  computer  engineering  at 
the  University  of  South  Alabama.  His  research  interests  include  smart 
energy  management,  ultrafast  computer  architectures  and  algorithms,  sig¬ 
nal/image  processing,  pattern  recognition,  and  digital  system  design.  He 
is  a  fellow  of  OSA,  a  fellow  of  the  SPIE,  a  senior  member  of  IEEE,  a 
member  of  ASEE  and  AIP 


